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Abstract. The interplay between the quantum interferences responsible for one particle localization over 
a length Li, and the partial dephasing induced by a local interaction of strength U with another particle 
leading to partial derealization over a length L2 > Li, is illustrated by a study of the motion of two 
particles put close to each other at the initial time. Localization is reached in two steps. First, before the 
time t\ necessary to propagate over Li, the interaction slows down the ballistic motion. On the contrary, 
after ti the interaction favors a very slow derealization, characterized by a log(f) spreading of the center of 
mass, until L2 is reached. This slow motion is related to the absence of quantum chaos in this one dimen- 
sional model, the interaction being only able to induce weaker chaos with critical spectral statistics. Under 
appropriate initial conditions, the motion remains invariant under the duality transformation mapping the 
behavior at small U onto the behavior at large U. 

PACS. 71.10.-w Theories and models of many electron systems - 71.30.+h Metal-insulator transitions 
and other electronic transitions - 73.20.Jc Derealization processes 



1 Introduction 

' In this third work of a series [1,2,3] concerning two inter- 
acting particles (TIP) in a one dimensional random lat- 
tice, we use the time dependent Schrodinger equation to 
describe the competition between the one particle quan- 
tum interferences induced by the random potential lead- 
ing to localization, and the mutual dephasing induced by 
a local interaction and leading to partial derealization. 

■ At the initial time, a wave packet representing two parti- 
cles in two neighboring sites is constructed in the middle 
of a disordered chain of size L. A repulsive on site in- 
teraction of strength U is considered. The particles are 
assumed to be two electrons with opposite spins, the ini- 
tial wave function is symmetric and remains symmetric 
during the quantum motion. We study the short times 
where the particles visit scales small compared to the one 
particle localization length L\ till the long times where 
the spreading of the center of mass saturates at the TIP 
localization length L 2 2> L\. For this, we use an efficient 
automaton-like algorithm adapted to discrete scalar wave 
propagation in a system of finite size L. We restrict the 
study to strong localization, such that the quantum mo- 
tion is the same when the size increases from L = 512 
to L = 1024 and t — > 00. This guarantees that our con- 
clusions are not biased by finite size effects coming from 
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successive boundary reflexions. The price to pay for this 
is to consider relatively small values for L\ and L 2 . 

The TIP-dynamics is characterized by two times t\ 
and t 2 , where a scale of order of L\ and L 2 is respectively 
explored. Between t\ and t 2 , we find that the spreading 
of the center of mass is extremely slow. To give an or- 
der of magnitude, L\ ss 16 can be quickly reached after 
t\ 100 (in units of time) while a time t 2 ~ 10 4 is nec- 
essary to reach only L 2 ~ 2L\. In this regime of interac- 
tion assisted propagation, we find that the center of mass 
spreads with a log(i)-law, quite different from a previously 
assumed diffusion law. This is consistent with the obser- 
vation [2] that the interaction can never drive the TIP sys- 
tem in one dimension to full quantum chaos with Wigner- 
Dyson spectral statistics. Only a weak critical chaos can 
be established, where the spectral fluctuations are statis- 
tically similar to those characterizing [4] different critical 
one particle spectra (3d Anderson model at the mobility 
edge [5] or certain pseudo-integrable billiards). Another 
important phenomenon illustrated by this study is the in- 
version of the effect of the interaction when the ballistic 
motion (t < ti) becomes sub-diffusive (t\ < t < t 2 ): first 
U defavors the ballistic propagation before having an op- 
posite delocalizing effect. 

The paper is organized as follows. The model and the 
used algorithm are introduced (Section 2). Then a series 
of useful results for understanding the complex features of 
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TIP dynamics (Section 3) are shortly summarized. Illus- 
trations of the TIP derealization phenomenon are given 
when t — ► oo. The role of the finite size effects are studied 
for L = 512, and the values of L\ and L 2 where they can be 
neglected arc estimated (Section 4). In the remaining part, 
we study how is reached this long time limit, after two suc- 
cessive regimes of the quantum motion. First, a ballistic 
regime (t < ti) is studied where U defavors propagation 
(Section 5). The dependence of the TIP dynamics on the 
chosen initial wave packet is illustrated. We consider both 
two particles put at t = on the same site (energy w U) 
and on two neighboring sites (energy 0). Following the 
initial condition, the dynamics probes two different sets of 
states in the large {/-limit: molecular states of energy w U 
and hard core boson states of energy s=a 0, as discussed in 
Rcf. [2] . The duality transformation mapping the behavior 
at small U onto the behavior at large U is illustrated for 
appropriate initial conditions. A study of t\ follows (Sec- 
tion 6) before describing the sub-diffusive regime where 
interaction favors a slow TIP propagation (ti < t < t 2 : 
Section 7). The log(t) spreading of the center of mass is 
related to the long life-time observed in Rcf. [1] for the 
free boson states (TIP states for U = 0). The relation 
between the observed very slow derealization, the multi- 
fractal measure [1] characterizing the interaction induced 
hopping terms coupling the free boson states, and the crit- 
ical weak chaos observed in Ref. [2] are discussed. 



2 Model and numerical algorithm 

To study the motion of two electrons with opposite spins 
in a one dimensional Anderson tight binding model with 
on site interaction, we have to numerically solve the dis- 
cretized TIP Schrodinger equation: 

i. mt +e))-m-e)))=nm)), a) 

where H = H ®1 + 1®H + H lat . H is the one particle 
tight binding Anderson Hamiltonian: 

L-l L 

#o = ]T (\n)(n + l\ + \n+l)(n\) + J2v n \n)(n\. 

n—l n—1 

L is the system size and V n are independent random vari- 
ables, uniformly distributed inside the interval [— W, W]. 
The ket \n) stands for the electronic orbital located at the 
site n of the one dimensional lattice. The eigenstates |a) 
of H are localized on a length L\ (with L\ w 24/ W 2 at 
the band center). H mt is the on site interaction: 

L 

Him = U^2\n)® \n)(n\ <g> (n|. 

n=l 

In what follows, all the lengths are given in lattice spacing 
units and all the energies in units of 1/e. The times are 
then expressed in the corresponding units (e). The sites 
are labelled from — L/2 to +L/2, such that the site n = 



is located in the middle of the chain. The initial condition 
corresponds to 

tp ni ,n 2 {0) = lp ni ,n 2 ( e ) = ;-y=(<*ni,o£n 2 ,po +^2,0^1,^0)1 

where (nin 2 \ip{t)) = ip ni ,n 2 (*)■ 

When e is small enough, the discrete time Eq.(l) has 
the same physical content as its continuous version. 

To solve Eq.(l), we use an automaton-like algo- 
rithm [6] , which relies on a formulation of discrete scalar 
wave propagation in an arbitrary inhomogeneous medium 
by the use of elementary processes obeying a discrete Huy- 
gens' principle and satisfying fundamental symmetries, as 
described in Ref. [6]. Our algorithm avoids the direct dis- 
cretisation procedure and incorporates the symmetries un- 
derlying the Anderson model at the lowest stage of the 
construction. As a consequence the algorithm preserves 
the unitarity of the dynamics, insuring the normalization 
of the wavefunction at all times, ^ „ 2 |VVu,n 2 (i)| 2 = 1> 
up to a small correction of order e 2 . Besides, the construc- 
tion is optimized for implementing the algorithm on mas- 
sively parallel machines. The numerical simulations have 
been carried for a time step e = 0.05. The simulations were 
performed on a 16K processor Connexion Machine. Given 
a value of the disorder strength W and a disorder config- 
uration, the wavefunction has been calculated for chains 
of length as large as L = 1024 and up to a maximum of 
10 6 units of time. 



3 Review of some useful results 

Since the two body quantum motion is quite complex, it 
is useful to have in mind a few previous results that we 
shortly summarize. 



3.1 Non linear a model 

The first analytical descriptions [7,8] of the TIP system 
are mainly based on simplified random matrix Hamiltoni- 
ans with independent Gaussian entries. The purpose was 
first to explain the new phenomenon of pair propagation 
at scales larger than L\ and the pair localization at a 
scale L 2 ^> L\. In Ref. [8], from such an effective random 
matrix Hamiltonian, a supcrsymmctric nonlinear a model 
was derived, closely related with the one found by Efetov 
for non interacting electrons in disordered metals. The ap- 
proach was mainly developed for an arbitrary dimension 
d. We recall the conclusions for a strictly one dimensional 
model (d = 1). 

Let us denote by \a) a one particle state located near 
the site n a , \a[3) the symmetrized product states forming 
the eigenbasis of the TIP Hamiltonian when [7 = 0. Since 
they are the symmetric eigenstates without interaction, 
we call them "free boson states" as in Ref. [2,3]. We de- 
note by r a p their inverse lifetime when the interaction is 
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switched on. The TIP supersymmetric a model gives A/3 
in agreement with the Fermi golden rule: 

A/3 ~ 2^|(a/3|A„t|a/3)| 2 z, off « A/ (^^) , ( 2 ) 

where 



B is the kinetic energy (band width), (a(3\H int \-fS) the 
typical interaction matrix element between two free boson 
states and u e s the density of free boson states |7<$) coupled 
to | a/3) by the interaction. /(x) w 1 when x < 1, i.e. when 
the typical distance \n a — np\ between the localized states 
| a) and \(3) is smaller than A and /(x) w exp(— x) when 
\n-a — np\ > A- The lifetime of |a/3) depends on the spa- 
tial overlap between |a) and ,3). Localized at a distance 
|«a ~ np\ large compared to A, the two particles have an 
exponentially small probability to be on the same site and 
to feel the on site interaction. The corresponding lifetime 
is exponentially large. For |a) and |/3) localized inside the 
same localization domain of size A, the lifetime A~ de- 
fines an important characteristic time. In Ref. [8], AT an d 
A are respectively the smallest resolved time and length 
scales. The main assumption of the approach, contained in 
the effective random matrix Hamiltonian from which the a 
model is derived, is that the one particle states are essen- 
tially ergodic, chaotic and random inside their localization 
domain. This is a simplification of the one body dynamics 
inside A, which amounts to assume for the single parti- 
cle a "zero-dimensional" dynamics inside A- For L = A, 
one gets (a{3\H int \jS) w ±U/L^ 2 and v e fi w L\/B near 
the band center, and eventually the above expression for 
A- As far as the dynamics are concerned, the main con- 
clusions of Ref. [8] are as follows. For t > A i t ne time 
where the pair size p{t) is roughly equal to \n a — np\ can 
be estimated from the lifetime of the free boson state |a/3) 
with \n a — np\ w p(t). This gives 

p(t)<x A(l + log(Ai))- (4) 

At the same time, the center of mass R(t) exhibit a dif- 
fusive motion R(t) = ^/D 2 (t)t, with a slightly time de- 
pendent diffusion constant D 2 (t) ~ (U 2 /B)(Li/ log(Ai))- 
This small time dependence of D 2 (t) comes from the fact 
that the frequency of the collisions between the two parti- 
cles decreases as the pair size grows. This diffusion stops 
when R(t) w A « {U/B) 2 L\, where TIP localization oc- 
curs. 



3.2 Level curvature 

Useful, though indirect information for the TIP propaga- 
tion at scales L < A can be found in Ref. [9] where the 
sensitivity of the TIP levels At to a change of boundary 
conditions is given. Detailed numerical calculations con- 
firming the predictions of Ref. [9] are given in the last 



paper of this series [3]. For a ring threaded by an AB- 

flux A the TIP curvature C 2 {E) = J2 A ^§^H e ~ e a) is 
given by the expression: 

C 2 (E)* gi ^--( gi -l)g 2 (U)^-. (5) 

gi and A\ are respectively the one particle conductance 
and mean level spacing. g 2 (L,U) can be understood as 
an interaction-assisted TIP conductance [10] of order 
r a p(L,U)/A2(L), where A 2 (L) is the spacing of the free 
boson levels directly coupled by the interaction. The above 
expression implies that the effect of U strongly depends 
on L. When L -C A> 9i ~ 1 ~ 9i, and C 2 is mainly given 
by a kinetic one particle term g\A\l A 2 reduced by a small 
correction proportional to g 2 (L,U). This means that for 
L <C A (ballistic one particle regime) the easy propa- 
gation due to kinetic terms is not yet strongly affected 
by the one particle quantum interference. In this case, 
the presence of the second particle defavors the propaga- 
tion of the first and the interaction slightly reduces C2, 
as confirmed in Ref. [3]. When L > L\, the one particle 
transport is suppressed by the quantum interferences (An- 
derson localization) and the term proportional to U in C 2 
changes its sign, since g\ — 1 w —1 up to exponentially 
small corrections. Thus, TIP transport is favored by the 
interaction, the presence of the second particle leading to a 
decoherence of the localizing quantum interferences of the 
first. C 2 is of order g 2 for the few TIP states re-organized 
by the interaction when L 3> A- From the behavior of 
the TIP curvature, one should expect a change of the role 
of U for the TIP quantum motion as a function of time: 
starting from a localized wave packet, the exploration on 
a scale L < A for t < t\ should require a longer time 
with interaction than without. But without interaction, 
the exploration of scales L 3> A is forbidden by Ander- 
son localization, while it becomes possible in the presence 
of interaction. An important issue is then to know what is 
the time scale t 2 — t\ required for this exploration, before 
the pair itself gets localized when t s=s t 2 - According to 
Ref. [8], t 2 — t\ = \J L 2 j D 2 {t 2 — t\). However, one should 
have in mind some further results obtained for strictly on 
site interaction and strictly one dimension. 

3.3 Multifractality 

The interaction matrix elements 

L 

(a/3|Ant|7<5) = 2C/^^*(n)^(n)A(n)AA), 

71=1 

(with A(^) = ("| a) ) defines a measure of the free boson 
states Y)$) coupled to a given |a/3). It was shown in Ref. [1] 
that this measure is multifractal. In contrast to earlier as- 
sumptions, the density v c s(L) of free boson states \j5) 
effectively coupled to |a/3) is much weaker than the total 
density v 2 (V) of free boson states. This is not surprizing 
when there is no disorder: only free boson states of same 
total momentum are directly coupled, and i/ c g — v\ oc L 
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and not v 2 oc L 2 . When disorder is switched on, the mo- 
mentum is no longer a good quantum number, and one 
gets v\ < v e g < v 2 . More precisely, if one wants to es- 
timate r a p using the Fermi golden rule, one needs the 
effective density of states |7<5) coupled by the second mo- 
ment (q = 2) of the interaction matrix elements. The mea- 
sure in the TIP Hilbert space of the support of this set of 
states is neither d — 1 (as for the clean case) nor d = 2 
(as for the chaotic d = case), but a fractal dimension 
1 < f(a(q = 2)) w 1.75 < 2. This gives a reduced effec- 
tive density v e s oc L/( a (9= 2 )) . The measure is multifractal, 
since this density depends on the considered q th moment of 
the coupling term. A direct implication of this multifrac- 
tal character is that the lifetime of the free boson states 
is longer than i~i for L as L\. One should expect that 
the diffusion law obtained ignoring multifractality under- 
estimates the time t 2 — t\ for the particles to propagate 
between L\ and L 2 . The multifractality of the set of di- 
rectly coupled free boson states should mean a very slow 
interaction induced derealization. 



3.4 Weak critical chaos 



The TIP spectral fluctuations also lead us to expect very 
slow dynamics before localization. For a given L, the TIP 
spectrum has Poisson statistics in the limits where either 
the disorder W or the interaction U are too weak or too 
strong. In the limit of the clean system (W = 0) , this is 
due to the fact that the effective density of coupled free 
boson states is v\ cx L <C v 2 oc L 2 . For W 3> 1, one has 
Li <C L and the small part of TIP levels being reorganized 
by U is totally hidden behind the main part of the non 
reorganized spectrum, corresponding to free boson states 
with \n a — np\ > L\ and which remain eigenstates when 
U is switched. For [/ « 1, the TIP states are basically 
the free boson states of energy e a + ep, i.e. almost un- 
corrected TIP-levels. For U » 1, the levels have again 
an energy e a + ep (duality property explained in the next 
paragraph). In the plane (U,W), inside those mentioned 
Poisson limits, in a domain centered around L w L\ and 
U w 1, the TIP spectrum becomes more rigid, though less 
rigid than a Wigner-Dyson spectrum associated to quan- 
tum chaos. The spectral rigidity saturates [2] to an inter- 
mediate rigidity between Poisson and Wigner-Dyson. This 
rigidity is however not arbitrary, but has a universal char- 
acter shared by many one particle 'critical' systems [4], 
such as the Anderson model at the mobility edge, a mixed 
system where integrable and chaotic trajectories coexist 
or a pseudo-integrable billiard where all trajectories be- 
long to a surface of genus larger than one. By the term 
"weak critical chaos" , we mean that under certain circum- 
stances (U w 1 and L a L\), the TIP system belongs to 
the same critical universality class than those one particle 
systems, at least as far as the spectral statistics are con- 
cerned. The interaction can never drive the TIP system 
towards a stronger chaos, i.e, towards full quantum chaos 
with Wigner-Dyson statistics. 



3.5 Duality 

The Hamiltonian without interaction is diagonal in the 
basis of the free boson states. As pointed out by Pono- 
marev and Silvestrov [11], there is an cigenbasis appro- 
priate when JJ — > co, having the same energies e a + ep 
than the free boson states. Indeed, when U — > oo, one has 
just to solve a non interacting problem with new bound- 
ary conditions. Since particles cannot be on the same site, 
we define L(L — l)/2 hard core boson states \hc) of com- 
ponents {n\n 2 \hc) given by 

^= [^aM^(ni) - **(ni)*p(n 2 )] ^~ n \ . (6) 
V2 \n 2 — ni| 

A hard core boson state is just a 2 x 2 antisymmetric Slater 
determinant resymmctrized by the factor (n 2 — nf) /\n 2 — 
n\\. To complete this basis of the symmetric TIP Hilbert 
space spanned by L(L + l)/2 states, we add L molecular 
states \nn) of energy 2V n + U, V n being the random po- 
tential of the n th site. This set of molecular states forms 
a small sub-band which goes to very large energies when 
U — ► oo . For the main sub-band of hard core bosons states, 
centered around E = as the free boson states, there is a 
duality relation [2,11] U <-» A/U, mapping the behavior at 
small U onto the behavior at large U . A = \p2A for E = 0. 
One finds that the coupling terms between the free boson 
states are given by 

n 

while the coupling terms between the hard core boson 
states arc a sum of terms like 



E 

n,n' ,n' 



^ a {n)W*{n)^{n')Mn") 
U + 2V n — E 



with various combinations of n', n" — n ± 1. The duality 
is obtained neglecting the difference between n and n' and 
assuming that U + 2V n — E m U. This gives a strong de- 
pendence of the TIP dynamics on the initial wave packet. 
For large U, the (\hc) U |nn)) form an cigenbasis, the mo- 
tion of a wave packet located at the sites and po for t = 
is given by the time dependent wave function: 



V'ni ,n 2 {t) 



E 

-(hc,nn) 



Cni,n 2 ,sCo,po,S e 



-iE s t 



where the summation S goes over the states \hc) and \nn) 
and C ni .„ 2j 5 = (ni,n 2 \S). For two particles put on the 
same site n at the beginning (i.e. with an energy U + 
2V n ), we mainly probe the molecular state \nn) when U 
is large: the two particles stay on the same site for a very 
large time. For smaller values of U, this initial wave packet 
starts to probe the molecular state \nn) before decaying 
onto the neighboring hard core bosons states of energy of 
order U. To eliminate the time motion associated with the 
molecular states and to see the time motion associated to 
the hard core boson states, which only exhibit the duality 
property, one has to start with two particles put close to 
each other on two different sites and po- 
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4 Asymptotic TIP localization and finite size 
effects 

The role of the repulsive interaction U on the quantum 
motion is shown in Fig. 1 and Fig. 2, illustrating the TIP 
derealization effect [7] in a strongly disordered chain for 
U = 1. This effect is a consequence of the mixing by the 
interaction of free boson states close in energy, delocalizing 
the TIP system in the free boson basis. Since the one body 
states are localized, this derealization in the free boson 
basis also means derealization in real space. The initial 
condition corresponds to po = 1- We have taken the site 
potentials Vo = Vi = 0, in order to benefit by the large 
density of TIP states for E w 0. 

In Fig. 1 and Fig. 2, the rainbow color code indicates 
on a linear scale the small values of \i>m,n 2 (t = 5. 10 4 ) | 2 in 
red up to the large values in violet. In the upper Fig. 1, 
one can see how the two particles are confined by the 
random potential without interaction, in a localization do- 
main which is very quickly reached (typically for t w 200). 
In the lower Fig. 1,17 = 1 and the center of mass becomes 
delocalized as sketched in the upper Fig. 2. TIP localiza- 
tion after ensemble averaging is shown on the lower Fig. 2 
for U = 1. This TIP ellipsoidal localization domain is 
reached and stops to spread after a considerably larger 
time (typically for t ps 5.10 4 ). For a given sample, one 
can see that \ip ni , n3 (t — 5.10 4 )| 2 does not homogeneously 
fill the ellipse, and is characterized by large fluctuations, 
mainly near the border of the ellipse. These fluctuations 
are somewhat similar to those characterizing the interac- 
tion matrix elements coupling the free boson states (see 
Fig. 1 of Ref. [1]). As we have checked, il) ni ,n 2 {t) develop 
an anisotropic multifractal behavior when the interaction 
assisted propagation begins to dominate the dynamics. 

To study the spreading R(t) of the center of mass and 
the size p(t) of the pair, we use the following functions: 



R(t) = 



256 




Pit) 



where m, 2 = J2 ni ,n 2 \^n y ,n 2 {t)?ni,2- We have checked 
that the disorder average of n 12 does not depend on time. 

Before going further, we first check that the localiza- 
tion effects are strong enough so that the TIP motion 
that we study corresponds to the dynamics of the infi- 
nite chain, and not of a finite chain with boundary effects. 
Indeed, for a too weak disorder, fast fronts of the wave 
function could propagate [12] up to the boundaries and 
be reflected. This will affect the long time behavior, as 
the reflected waves enhance the localization of the center 
of mass. The enhancement factor for R(t) would be under- 
estimated. Since our numerical algorithm allows to study 
systems of size up to L = 1024 for 10 6 steps of time, we 
compare in Fig. 3 the motion for L = 512 and L = 1024 
and for L\ — 16 and L\ = 36. 



-256 




-256 



256 



256 



n 2 



-256 




-256 



Fig. 1. |^m,ni 
and po = 1 at t 
(7 = 1. 



256 



(t = 5.10 )| for two particles put on sites 
= 0. L = 512 and Li = 16. Up: U = 0. Down: 



For Li — 16, R(t) saturates at the same value L 2 « 36 
when the size L = 512 is doubled, while R(t) has strong 
finite size effects for Li = 36 above t 10 3 when L = 512. 
The study of TIP derealization in the infinite chain can 
be investigated when L = 512 only for a large disorder, 
where both L\ and L 2 are relatively small. Since L 2 is 
reached after a very long time, our numerical method is 
not convenient to study how L 2 depends on L\ and U. 
We just show on Fig. 4 the probability density p(R) = 
J2 ni+ n 2 =R \Tp ? m 2 (t = 5.10 4 )| 2 for U = and U = 1 when 
L\ = 16. This proves that the center of mass is indeed 
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256 



■256 




■256 



256 



Fig. 3. Finite size effects for (7 = 1. R(t) for L\ — 16 and 36 
for two sizes L — 512 and 1024. 
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Fig. 2. Up: Scheme of the TIP derealization effect: 
\4>ni,n 2 (°o)\ 2 is concentrated in an ellipse corresponding to a 
center of mass R delocalized by the interaction on a scale larger 

than the pair size p. Down: same as in Fig. 1 for U — 1, L 1 = 36 Fig. 4. TIP derealization effect. Probability distribution 



after having averaged over 20 samples 



p(R(t = 5.10 4 )) for a single sample with Li = 16, L = 1024, 
U = (dashed line) and U — 1 (thick line). 



exponentially localized over a length Li ss 2Li, without 
finite size effect. For larger values of L\ we have seen larger 
effects (R(U = 1)/R(U = 0) « 3.5) at long times, but this 
only gives a lower estimate for the enhancement factor 
L2/L1, boundary effects being non negligible. 

The different regimes of the quantum motion 

We now study the intermediary time scales during which 
the center of mass R(t) spreads, before the time where 
it saturates and TIP localization occurs. For U — 0, the 
aspect ratio of |^m,n 2 (i)| 2 ) defined by R(t)/p(t), remains 
equal to one at all times, but for U =/= 0, the time evolution 



of this ratio exhibits three regimes (Fig. 5), delimited by 
two characteristic time scales t\ and t%. For t < t%, the re- 
pulsive interaction favors p(t) and defavors R(t). The ratio 
R(t)/p(t) decreases. This is the ballistic regime character- 
izing the length scales smaller than L\. The situation is 
opposite for t\ < t < t% where L\ has been reached and 
the interaction assisted propagation of the center of mass 
begins, on scales larger than L%. The increase of R(t) is 
now much faster than the increase of p(t), and the ratio 
R(t)/p(t) increases. Li is reached at t = t% where TIP 
localization occurs. 
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Fig. 5. Single sample with L\ = 16, L = 1024 and (7 = 1. Up: 
R(t)/p(t). Down: R(t), p(t) and p(t) for U = and [7 = 1. 

5 Ballistic propagation and duality 

For t < t\ we find that the spreading of the center of mass 
is almost ballistic: 

R(t) ~ v{U)t^ u) with n(U) w 1, 

and that the interaction reduces the increase of R(t). The 
time evolution strongly depends on the initial condition. 
When the two particles are injected on the same site at t = 
0, with an energy of order U, the spreading of the center of 
mass is almost suppressed by a too large interaction. This 
is the dynamics associated to the molecular states \nn), 
which do not decay when U becomes very large (Fig. 6). 

On the contrary injecting the two particles at two 
neighboring sites (0 and po = 2), one can see the dy- 
namics associated to the hard core boson states and the 
consequence of the duality relation U <-> l/U between the 
free bosons and the hard core bosons. Since the density 
of those states exhibits a Van Hove singularity at E = 



Fig. 6. R(t) at small times for different values of the inter- 
action. Li = L = 200, average over 20 samples. Up: po = 0. 
Down: po = 2. 



for [7, W = 0, we have also taken Vq = V Pa = 0. This op- 
timizes the coupling between ip(t — 0) and the hard core 
spectrum. The duality shows up in the quantum motion. 
On Fig. 6, one can sec (i) that when p(t = 0) = 2, R(t) 
increases almost linearly as a function of t (i.e. that the 
motion is almost ballistic, (ii) that U < 2 decreases R(t), 
(iii) that R(t) is very similar when U = and [7=15 
(duality) . 

On Fig. 7, the averages over the random potential of 
R(t = 40) and R(t = 1000) are shown as a function of 
the strength U of the interaction. The upper figure corre- 
sponds to the ballistic regime where U defavors transport, 
while the lower figure corresponds to the case where U has 
started to favor transport. The two curves illustrate the 
consequences of the duality between the free bosons and 
the hard core bosons and the inversion of the role of U 
around t\. 
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u 



Fig. 7. Duality and inversion of the role of U: R(t) averaged 
over 100 samples. Up: ballistic regime t = 40 < ti, L = Li = 
200. Down: sub-diffusive regime, t = 1000, L = 512, Li = 36. 



6 Scaling in the ballistic regime 

To define the characteristic time t\ separating the ballistic 
and sub-diffusive regimes, one can use many criteria. For 
instance, t\ can be defined (i) as the lifetime h/Ti of the 
free (hard core) boson states, when U < U c ~ 2 (U > 
U c ): (ii) as the time if where the particles reach the one 
particle localisation length L\ when R(U, tf ) = R(U = 
0,t = oo); (hi) as the time i™ 1 " where the minimum of 
R(t)/p(t) is reached (see Fig. 5); (iv) as the time scale if 
allowing to map the curves R(t)/p(t) onto a single scaling 
curve R(t)/p(t) = f s (t/t{). The existence of such a scaling 
and the L\ and U dependence of if are shown on Fig. 8. 

We have checked that the definitions (iii) and (iv) are 
compatible. For t < t\ the motion is essentially ballistic, 
and we expect that t\(L\) oc L\, as for a clean system of 
size L — L\ (definition (i)) where a term of order ±U/L\ 



30 




12 3 

U 



140 
120 




t/t, 1 



Fig. 8. Up: t? in (U) with L x = 2A,L = 256 and 200 samples. 
The line is a linear fit: ti = Li - SU. Middle: C in (Li) with 
L fa 3 x Li, U = 1 and 50 samples. The line is a linear fit: 
tf ln = 0.5Li. Down: rescaling R{t/t{) with Li = 24, U = 0.5 
(pluses), [7 = 1 (squares), U = 1.5 (diamonds), U = 2 (full 
triangles) and U = 1 Li = 16 (triangles), Li = 36 (stars), 
Li = 50 (crosses). 

coupled v\ oc Li free boson states. 

= f(U) w 1 - 0.3(7. 

The interesting feature of t™ ln is that i? w p at this time, 
as when [7 = 0. This time where there is an inversion of 
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the effect of the interaction should be related to the size L 
where the TIP level curvature [3] does not depend on U. 
According to Rcf . [3] , this size is of order (but not exactly) 
Li. 



7 Very slow derealization and weak critical 
chaos 

After the ballistic propagation for t < t\, the spreading 
of the center of mass measured by R(t) saturates without 
interaction. This is due to one particle quantum interfer- 
ences yielding one particle Anderson localization. When 
U 7^ 0, this saturation is suppressed, but the spreading 
R(t) has now a so slow increase that a logarithmic scale 
for the time t is appropriate. 

Let us first consider the increase of the relative sep- 
aration p(t) between the two particles when interaction 
assisted propagation begins. As recalled in the subsec- 
tion 3.1, we expect a behavior given by 

p(t)«Li(l + KTii)), 

which turns out to be in good agreement with the nu- 
merical results. Plotting p{x) — p(ti) as a function of 
x = log(t/ii), one can check on Fig. 9 the predicted loga- 
rithmic behavior. 



100 



10 





// 
/ / i* 

/ / / 




/ '/ / 




L,=16 






L=24 






— L,=36 






— y~x 


/ // ■'' 




/ // / 



0.01 



0.10 1.00 
x=Log(t/t,) 



10.00 



Fig. 9. Dynamics of the size of the pair: p(x)— p(t\) as a func- 
tion of x = \og(t/ti) in log-log coordinates for L\ = 16, 24, 39, 
L = 512 and U = 1. 



The evolution of the center of mass is on the contrary 
not described by the (modified) diffusion law 



R(t) « y/D 2 (t)t, 

but has a much slower motion. However, this slow inter- 
action induced delocalization is not very surprising since: 
(i) The hopping terms induced by the interaction be- 
tween the free boson states (or the hard core boson states) 



are much smaller in one dimension than assumed in the 
random matrix model used to derive the non linear a 
model of Ref. [8] (subsection 3.1). The interaction ma- 
trix is not Gaussian and the effective density of coupled 
free boson states is significantly reduced (subsection 3.3). 
Due to this multifractal perturbation, the lifetime of the 
free boson states is much larger than for a more normal 
(Gaussian) interaction matrix, as checked in Ref. [1]. 

(ii) For L w L\, the interaction can only drive weak 
critical chaos (subsection 3.4). For L > L\, the spec- 
trum becomes less rigid, with statistics intermediate be- 
tween critical statistics and Poisson statistics. Indeed, for 
L > L\, the TIP spectrum becomes a superposition of 
many states not reorganized by U and having uncorre- 
cted fluctuations (Poisson) and a small part of states 
having critical statistics (weak critical chaos). With the 
chosen initial condition (two particles put close the one to 
the other at t = 0) one can argue that we mainly probe 
the few states with critical spectral statistics. Those criti- 
cal statistics are associated with slow anomalous diffusion. 
Let us take one of those billiards with critical statistics: 
a right triangle [4] with smallest angle equal to ir/5 and 
Dirichlct boundary conditions. For a very long time, the 
classical trajectories are stable, the system remains on the 
same KAM torus in phase space, until the corner with an- 
gle 47r/5 is reached. At this moment only, the trajectories 
can escape from the original KAM torus and start to ex- 
plore other parts of the phase space. This suggests us a 
possible analogy between a single particle in the triangle 
billiard and the two particles in a disordered chain of size 
w L\ with on site interaction: each particle is trapped in 
the one particle phase space between the collisions. Only 
after a collision, the frequency of those collisions being 
very low and depending on the strength of the random 
potential, the two particle phase space starts to be ex- 
plored. This may explain the similarity with a single par- 
ticle in the triangle. Critical statistics are associated with 
very slow explorations of the phase space, and hence of the 
real space. The fact that the interaction can never drive 
full quantum chaos, but only weak critical chaos, makes 
likely a very slow interaction assisted propagation. If the 
motion was isotropic in the plane (711,712), we should ex- 
pect anomalous diffusion (R 2 (t) w t a ) with a < 1, as 
observed for a single particle at the mobility edge in three 
dimensions, where the spectral statistics are critical. 

Since the motion is anisotropic in the plane (711,772), 
we do not find simple anomalous diffusion, but a simple 
log(£) behavior: 

R(t) cx log(i), (7) 

as shown on Fig. 10 for L\ — 16 (see also Fig. 3 for L\ = 
36). 

This logarithmic delocalization is not easy to explain, 
but reminds us the classical problem of a random walk 
in a percolating network. Without anisotropy, one gets 
anomalous diffusion. When one introduces anisotropy in 
the walk, it has been observed [13] that the dynamics is no 
longer described by a power law R(t) w t a , but by a log(i) 
law. If we order the localized one particle states |a) by the 
location n a of the center of their localization domain, from 
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Fig. 10. 

U = 1 



Slow derealization of R(t): L\ = 16, L = 1024 and 



one side of the chain to the other, the TIP model in the free 
boson basis \a(3) is both anisotropic (larger hopping along 
the center of mass direction than along the other direction) 
and the hopping terms (a(3\Hi nt \jS) have a multifractal 
measure. One can then argue that the two problems might 
be related, sharing the same log(t) spreading of R(t). 



tian Vanneste and Dietmar Weinmann for their con- 
stant interest in this work. 
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8 Conclusion 

In summary, this study of the TIP quantum motion has 
given some new insights on several aspects of the model. 
There are two length scales L\ and L 2l and two corre- 
sponding time scales t\ and ti. For L < L\ (t < ii) the 
interaction defavors the pair propagation and reduces the 
level curvature. On the other hand, for L > L\ (t > t\), 
there is a very slow interaction assisted pair propagation 
and the level curvature increases as a function of U. Our 
main result is that this derealization is very slow, in qual- 
itative agreement with the concepts of interaction induced 
weak critical chaos and of multifractal hopping terms. 
Moreover, the wavefunction \^ ni ,n 2 (t)\ itself has multi- 
fractal features visible in Fig. 1. With appropriate initial 
conditions, one can observe the symmetry U <-» A/U. In 
conclusion, we underline that our results characterize sym- 
metric states with purely on site interaction in strictly one 
dimension. It will be interesting to study if they remain 
valid for longer range interactions, or in a quasi-one di- 
mensional limit, or if the behavior predicted from the a 
model approach of Ref. [8] becomes valid. 
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